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Distillation protocols enable generation of high quality entanglement even in the pres- 
ence of noise. Existing protocols ignore the presence of local information in mixed states 
produced from some noise sources such as photon loss, amplitude damping or ther- 
malization. We propose new protocols that exploit local information in mixed states. 
Our protocols converge to higher fidelities in fewer rounds, and when local informa- 
tion is significant one of our protocols consistently improves yields by 10 fold or more. 
We demonstrate that our protocols can be compacted into an entanglement-pumping 
scheme, allowing quantum computation in distributed systems with a few qubits per 
location. 

1. Introduction 

Maximally entangled pairs of qubits are a powerful resource for a plethora of com- 
putational tasks that rely on quantum correlations, from quantum cryptography to 
quantum computation. Though useful, entanglement is a delicate property that is 
difficult to produce and store. Imperfect experimental conditions and noise from 
the environment produces statistical ensembles of quantum states that possess 
non-maximal entanglement. Fortunately for our aspirations to develop quantum 
technologies, entanglement can be distilled from many imperfect base pairs into a 
smaller number of more entangled pairsQ 

Whilst the literature on entanglement distillation is vast, it almost exclusively 
tailors to distillation of mixed states that are ensembles of maximally entangled 
statesJUMHH Such mixed states are said to be locally maximally mixed, as local 
measurements on one qubit have completely random outcomes. For all other mixed 
states there exist local measurements with non-zero expectation value, and so we 
know some local information about the mixed state. Some protocols propose to 
distill base pairs by destroying this local information in a process called twirlinj^, 
and other protocols simply ignore the local information^. But surely knowledge is 
power and no good can come from consigning local information to the flames. This 
intuition is the inspiration for this study of how we can exploit local information 
to design more streamlined distillation protocols. 

^Throughout we only consider two qubit mixed states. As is common in the literature, we use 
pair to refer to two qubit states. We reserve the phrase base pair for the entanglement resource 
obtained before any distillation has been performed. 
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After describing our protocols, we consider noise models grounded in experi- 
mental considerations, and compare the performance of various protocols at dis- 
tilling these base pairs. Any qubit embodied in the non-degenerate energy levels 
of a physical system will undergo amplitude damping and thermalization with the 
environment, which leads to a preference for the ground state and gives us local 
information. Also, the process of producing an entangled base pair can inherently 
produce local information. For example, the effect of photon loss in certain optical 
based procedures for entanglement generatiorP, where the mixed state contains an 
incoherent component of the doubly emitting state |1, 1). We will analyze different 
approaches to distilling (a) combined photon loss and depolarizing noise; and (b) 
amplitude damping noise for a range of different relaxation bases. 

Our protocols are best understood as building blocks. Having applied one round 
of a protocol, higher fidelities are achieved by repetition, which can be achieved in a 
variety of waydmHEHl e.g. by symmetrically always pairing up identical noisy Bell 
pairs, or by continuously pumping with the same base pair. The problem of iden- 
tifying the best approach to repetition can be viewed as a self-contained problem, 
which has been called the scheduling problenPSl We calculate the yield when per- 
forming distillation between identical pairs, so-called symmetric scheduling. Here 
we find that one of our protocols consistently attains higher yields — roughly 10 
fold or greater — in the strong local information regime. 

Since we are primarily concerned with improvements on the level of the building 
blocks, we will not give detailed comparisons of different approaches to scheduling. 
However, it is important that our protocols can be adapted for entanglement pump- 
ing schedules that conserve spatial resources at the cost of temporal resources. In 
section [5] we extend our protocols to enable entanglement pumping, where the pairs 
are non-identical. We will discuss the application of entanglement pumping to dis- 
tributed quantum computers with a few qubits at each local sit^331H Our pumped 
protocols demonstrate that all possible entangled base pairs can enable quantum 
computing with only two or three qubits per location, a result previously known 
for only locally maximally mixed states. 

Before embarking on a description of our protocols, it is first necessary to review 
some of the existing approaches to entanglement distillation, and also to introduce 
some technical tools for understanding the nature of local information and its rela- 
tionship to entanglement theory and local filtering theory. 



1.1. The prior art of entanglement distillation 

A density matrix is locally maximally mixed when all local observables have an 
expectation value of zero, otherwise the density matrix contains local information. 
It is well known that locally maximally mixed (LoMM) states can always be locally 
rotated so the density matrix is diagonal in the basis of Bell state^if. These states 
are entangled if and only if they have a maximally entangled fraction, -Fmax; greater 
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than 1 /2, where: 



max — max; 



(1) 



which is maximized over all maximally entangled | ■(/')• Bennett et al. proposed the 
first protocol for distilling these base pairs ^ ^, which was later improved by Deustch 
et al.^. The improved protocol is often called DEJMPS after the initials of its 
originators, and a single round proceeds as follows: 

(1) take two copies of p and locally rotate]^ such that the largest maximally entan- 
gled state is mapped to a particular Bell state, \^) |(f>+) = (|00) + |ll))/\/2; 

(2) perform local bilateral control-not gates, CX^^CX§^; 

(3) measure qubits Al and Bl in the computational basis; 

(4) we have succeeded and keep the final state if we measure |0)yii|0)_Bi or 
|1)ai|1)bi, otherwise we must restart with a fresh pair of p. 

For completeness, we define the familiar control- not gate: 



with X2 as control and xi as target. Throughout, we label Al and Bl as the qubits 
sharing one copy of p and A2 and B2 as sharing a second copy of p, with A and B 
denoting the location of the qubits. After a successful round, the fidelity, w.r.t l^"*"), 
and success probability of DEJMPS depend on only the diagonal elements of p in 
the Bell basis. If the initial fidelity exceeds 1/2, then a successful round produces 
a final state with an increased fidelity. To achieve higher fidelities DEJMPS must 
be repeated, with different approaches to repetition discussed in sections |4] and [5] 

DEJMPS can also be used to distill mixed states that are not of LoMM form, 
but its performance has no dependence on terms that are off-diagonal in the Bell 
basis. As such, DEJMPS does not utilize any local information that may be con- 
tained in these off-diagonal terms. When applying DEJMPS to non-LoMM states, 
it still requires -Fmax > 1/2 for the protocol improve fidelity towards unity. How- 
ever, there are many entangled states that this protocol fails to distill because 
Fmax < 1/2. In these instances, ignoring the off-diagonal terms effectively destroys 
the entanglement along with the local information. 

Though the literature has focused on LoMM states, a paper by the Horodecki 
family confronts the problem of non-LoMM states and proves that all entangled 
states can be distillecP^l. in short, the Horodeckis found a method of taking any 
entangled state and probabilistically producing a state with i<"niax > 1/2, hence 
enabling DEJMPS to be applicable. The Horodeckis use the idea of local filtering, 
where a local Positive Operator Valued Measurement (POVM) is performed and 

''This local rotation is not explicitly prescribed in Ref. |6] We have found that fidelity is greater 
when we diagonalize the 3 by 3 submatrix of the Ri j matrix for i,j = 1, 2, 3 (see Eq. Q). We also 
always apply local operations that interchange the Bell state components to maximize fidelity. 
These slight improvements on the original protocol have a modest effect on increasing fidelity and 
yield. 



(2) 
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we post-select on one outcome. Formally, we denote a local filter on qubit A by 
/a, which can be any one qubit matrix with singular values of magnitude less than 
onc|^ The post-filtering state is: 

^ if A ® tB)p{fA ® IbV f . 

^™ tr((/^®lB)p(/A®lB)t)' ^' 

where the denominator is the probability of successfully applying the local filter, 
which we denote P^^°'^- We will not focus on the communication cost of protocols, 
but note that every local filtering operation requires one classical bit of information 
to communicate success. The Horodeckis have an elegant proof that for all entangled 
states there exists a local filter, which we call the Horodecki filter, that boosts 
-Pmax above the 1/2 threshold. However, in the very same paper the Horodeckis 
identified a class of states for which the Horodecki filter works but is very inefficient. 
For this class of states they also identified a different filtering operation that was 
superior to the Horodecki filter. Consequently, we are left with the open question 
"is there a systematic and efficient procedure for exploiting local information to 
distill entanglement?" . 



1.2. Local information, local filters and concurrence 

Here we review some results on the interplay between local information, local fil- 
tering and entanglement. We begin by noting that local information of a density 
matrix is most apparent in its Hilbcrt-Schmidt decomposition: 

where (To_i^2,3 are the Pauli matrices 1,X,Y,Z and tr^p) = 1 entails _Ro,o — 1- 
If wc measure an observable <7i (8" cTj, then the expectation value of the measure- 
ment outcome is tT(p(ai (E) Cj)) — Ri_j. When measuring an observable local to 
qubit A, all the relevant information is encapsulated by Ra = (^f^i,Oi ^2,0i ^3,o)) 
and similarly local measurements on qubit B have expectation values predicted by 
Rb = (Ra.ii Ro.2, Ra.s)- Representing this information as vectors Ra and Rb is the 
Bloch sphere representation of each qubit. Hence, these six numbers describe all the 
local information, and when they all vanish the state has no local information and 
is a LoMM state. All other entries of the matrix, Rij {i ^ 0, j ^ 0), describe cor- 
relations between the two qubits. One advantage of the R matrix representation is 
the simplicity of how it transforms under local unitaries. Applying a unitary Ua to 
qubit A is equivalent to a rotation in the Bloch sphere, and describing this rotation 
by a proper orthogonal matrix Oa, it follows that Ra ~+ OaRa- Similarly, a local 
unitary Ub has an associated Ob, such that the unitary effects the transformation 



■^This must be true because otherwise there will exist density matrices for which the corresponding 
POVM measurement has a success probability of greater than 1. 
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Rb ObRb- When the joint system is rotated by Ua ' 
transforms as: 



U B , the whole R matrix 



R 



(5) 



1 \ /I 
QOa) \QOb, 

Verstraete et aZ.^ studied how local filters transform the R matrix, and found that 
local filters caused Lorentz boosts on R. They then showed that all physical states 
can be locally filtered into one of two canonical forms, one of which is LoMM: 

/I \ 

Ri.i 



MM 



Ri,i = +Pi 

R2,2 = -Pl 
R3,3 = +Pl 



P2 
P2 
P2 



P3 
P3 
P3 



Pa; 
Pa; 
Pa; 



(6) 



R 



R.D 



(7) 



i?2,2 
VO i?3,3. 

where pi are the probabilistic weights of 4 Bell states. The second canonical form 
they discovered they called the rank deficient (RD) class, as prd always has rank 
less than the maximal 4. The canonical form of RD states is: 

/I c \ 
a 
-a 
V6 l + h-cj 

where a, 6, c are reals. This form is not unique as the local information can be 
increased or suppressed by further local filtering. However, the local information 
can never be completely filtered away. Having already remarked that Horodecki 
filtering is not optimal, this motivates considering whether local filtering into the 
appropriate canonical form will perform better. Indeed, this will prove to be a good 
starting point. However, the class of rank deficient states is measure zero and so 
real world mixed states are unlikely to neatly fall into the rank deficient class. On 
the other hand, all full rank states can be local filtered into LoMM form, but if they 
contain a large amount of local information then the local filters will rarely succeed 
(see sections [3] and |4] for details). We will propose an approach that circumvents 
these problems. 

It is also instructive to consider how the entanglement of mixed states changes 
after local filtering. As an intermediate step to calculating the concurrence^, we 
find the eigenvalues of a matrix: 



where the tilde operation is 



p={Y ®Y)p*{Y ®Y), 



(8) 



(9) 



with the star indicating complex conjugation in the computational basis. This ma- 
trix will have 4 eigenvalues, labeled A,j and the concurrence is the largest eigenvalue 
minus the other 3. It is intriguing that local filtering may change the concurrence 
whilst preserving the ratios of concurrence eigenvalues. Formally, if we denote the 



Stepl PLoMM/Zinf Steps 2-4 

AlO-^—OBl > O O 1 ^ O-^i-O 

^20— r-052 ^ O O 

P PLoMM/Zinf 
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Fig. 1. A schematic of the LoMM and Zinf protocols. Both follow the same general pattern, but 
differ in the operation used in each step. They both consume base pairs, p, and output a Bell 
diagonal mixed state p2- 



post-filtering concurrence eigenvalues by A^, we find AJ/A^ = Xi/Xj. This property 
is useful as it entails that the pi of any post-filtering i?LoMM follow directly from 
the concurrence eigenvalues of the original p via: 

K = A,/(Ai + A4 + A3 + A2). (10) 



2. Two Proposed Protocols 

In this section we propose two protocols for entanglement distillation, called LoMM 
and Zinf protocol for reasons that will be explained. The LoMM protocol has never 
been explicitly proposed, but one could arrive at it by combining the single copy 
local filtering protocols of Rcfs. E] and [H with DEMPJS. Whilst LoMM will prove 
useful in some regimes, it can suffer from a pathologically small success probability 
when the initial mixed state contains a lot of local information. Our more elaborate 
Zinf protocol is quite unique compared to previous protocols. It is tailored to avert 
pathologies in the strong local information regime, and hence will outperform other 
known protocols in this regime. In this section we will describe one round of these 
protocols (outlined in Fig. fl]), with one round performance analyzed in section Is] 



2.1. LoMM protocol 

As indicated earlier, we begin by considering distillation via locally filtering mixed 
states into their canonical form, as defined by Verstraeteei al.-^. Typically, mixed 
states will local filter into LoMM form and so we call this the LoMM protocol. 
Given a base pair, p, we apply a local filter /a fs to bring the state into LoMM 
form, such that 

which succeeds with a probability, denoted -PlqMM' equal to the denominator of 
this expression. By conservation of ratios of concurrence eigenvalues under local 
filtering, we conclude this mixed state is: 



PLoMM 



= (Ai|$+)($+| + A2|^-)(*-| + A3|*+)(*+| + A4|$-)($-|)/^ A, 
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where are the eigenvalues of the concurrence matrix of p. Since local unitaries 
can interchange Bell states, we can choose Ai and A4 to be the largest and smallest 
Ai, respectively. After two successful attempts at local filtering, p — s- ploMM, we 
have two copies of ploMM- Now we have two copies of a Bell diagonal state it can 
be distilled in the usual fashion, and so the whole LoMM protocol is: 

(1) repeatedly attempt to local filter copies of p into LoMM form until we have two 
copies of PloMm; 

(2) perform local bilateral control-not gates, CX^-^CXgf; 

(3) measure qubits Al and Bl in the computational basis; 

(4) we have succeeded and keep the final state if we measure |0)yii|0)Bi or 
|1)ai|1)bii otherwise we must restart with a fresh pair of PloMM- 

The probability of success for steps 2-4 is: 

^^L^MM = ((Al + + (A3 + A2)^) / J2 (12) 

i 

and the final state is: 

P2 cx (A2 + a2)|$+)($+|+2AiA4|$-)($-| (13) 
+ (A2-hA^)|*+>(*+| + 2A3A2|*-)(*-|. 

Recall that we had the foresight to fix Ai and A4 as the largest and smallest Ai, so 
that the maximally entangled fraction: 

F' = ^1+^4 

'"^'^ (Ai+A4)2 + (A3 + A2)2' ^ ' 

cannot be increased by any reordering of the A^. Note also, that entanglement 
increases whenever the base pair was entangled. 

Whilst the LoMM protocol certainly has it merits, in later sections we will see 
that Plomm t)e extremely small. Before turning to this analysis we will first 
introduce our second, more sophisticated, protocol. 



2.2. Zinf protocol 



Since it is already known that we can always find local filters that bring the Hilbert- 
Schmidt decomposition into a form -RloMM or Rrd, it follows that all states can 
be locally filtered to some linear combination of these possibilities: 



/ 1 i?3,o\ 



i?2,2 

Vi?0,3 R3,3j 



(15) 



•^We reorder A; by applying the appropriate local unitaries to PloMM • 
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where the label "Zinf ' emphasizes that the only ofF-diagonal contributions of Rzmt 
are from the local Z-in/ormation. Our second protocol prescribes local filtering 
into Zinf form. However, for any given p there is a whole family of appropriate local 
filters that perform a transformation p — s- pzinf • Any local filters that belong to this 
family we denote gA® Qb, so: 

^ {9A ® 9b)p{9a ® gsV j-^g^ 

Recall that the denominator gives the local filter's success probability, PzlnT - If we 
already know that local filter Ja ® Jb brings p in LoMM form, then the family of 
gA 'S) gB local filters is: 

gA®gB = {ZA ® Zb){ua (X" UB){fA «> Jb), (17) 

where ua <8) ub are local unitaries that interchange Bell states, and za ^ zb are 
computational basis local filters such that: 

cxc,|0)(0| + (l-c,)|l)(l|, (18) 

with real ca and cb- Hence, this family of local filters has a simple description. 

Given this family of local filters, we are free to choose the local filter according to 
any figure of merit we desire, which will vary depending on the ultimate application 
of the distillation protocol. In the forthcoming analysis we select the local filter that 
maximizes the overall success probability: 

pall / pfiltcr\2 pdistil ^10^ 

-'Zinf — l-f^^Zinf ) Zinf • V^^J 

This is prudent as it can significantly boost the success probability above that 
attained by the LoMM protocol. Amazingly, we will see that whatever pzinf we 
choose to locally filter to, our distillation protocol will output the same final state 
up to an interchange of the concurrence eigenvalues. At this juncture it is timely 
to explicitly expand pzinf into bra-ket notation: 

pzinf - +P2|*")(*"| +r (|$+)($-| + |$-)($+|) (20) 

where pi follow from Eq. ^ and the off-diagonal terms relate to the local informa- 
tion, via r = {Rsfi + i?o,3)/4 and s — {Rs^o ~ Ro.s)/^:- For subsequent reference, we 
note that the concurrence eigenvalues of pzinf are: 



ipl+pl)/ 2-r^; 
{Pl -P2)yJ{Pl +P2)/2-r2; 
{pI+pI)/2 -s^- ^^'> 

(P3 - P4)V(P3 +P4)/2 - 

Once two copies of pzinf have been successfully prepared, we are ready to perform 
distillation. However, since pzinf is not LoMM we have to adapt the usual approach: 



Ai 


= VA + 


B; 


A ^ 


A2 


= Vc- 


D; 


B = 


A3 


= Vc + 


D; 


C = 


A4 


= VA- 


B; 


D = 
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(1) repeatedly attempt to local filter copies of p into Zinf form until we have two 
copies of pzinf; 

(2) perform local bilateral control-not gates, CX^^CX§f; 

(3) measure qubits Al and Bl in the computational basis; 

(4) we have succeeded and keep the final state if we measure bi, otherwise 
we must restart with a fresh pair of pzinf ■ 

The main difference to conventional distillation is that we post-select on only the 
|1)ai|1)_bi measurement outcome, which occurs with probability: 

i'ztf = l{iPi+P2r + {P3+Pi)' 4(r2 + s^)), (22) 

= ^((Ai + A4)2 + (A3 + A2)2). 

Notice that the success probability is upper bounded by 1/2, because we post- 
select only one measurement outcome. In contrast, both DEJMPS and LoMM post- 
select on two possible outcomes, e.g. |1)ai|1)bi or |0)yii|0)Bi, which allows them 
to approach unit success probability. Although this is a disadvantage for Zinf, this 
more restrictive post-selection can cope with the presence of local Z information and 
hence proves invaluable in distilling states with large amounts of local information. 
Upon success A2 and B2 are projected into the final state: 

P2 cx HpI + pI)/2 - r2)|$+) ($+1 + (pip2 - r^)\<^-){^-\ 
+ {{pI+pI)/2 - s^)\^+){^+\ + (P3P4 - 



which can be simplified using Eq.s (21 1 



P2 oc (A^-f a2)|$+)($+|+2AiA4|$")($-| (23) 

+ (A2 + a2)|*+)(*+| + 2A3A2|*-)(*~|. 

Again we see that Zinf produces a final state that appears to depend on only the 
concurrence eigenvalues, seeming to imply that Zinf always attains the same fideli- 
ties as LoMM. However, there is a subtle difference between the fidelities achieved 
by these protocols. For the LoMM protocol we applied local unitaries to ploMM 
ensuring that Ai > A2, A3 > A4. Here, the local unitaries, ua®ub, that interchange 
Bell states are buried in the middle of gA® 93, and consequentially they can affect 
the filter success probability. If we local filter to maximize -Pfj^f , then sometimes 
the Ai will be optimally ordered and sometimes not. Note that we cannot simply 
interchange Bell states after applying gA® 9b as applying local unitaries to pzinf 
can change local Z information into local X information. 



3. One Round Performance 

We have reviewed the DEJMPS protocol and Horodecki distillation and proposed 
our own protocols, LoMM and Zinf. We now turn to a comparison of the "one 
round" performance of these protocols for two families of physical noise model. A 
"many round" analysis will follow in subsequent sections. 
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Error rate: e Error rate :e 



Fig. 2. One round performance of distillation of photon loss noise (see Eq. l |24| ) for various 
protocols. Plots (al) and (a2) have the error rate e fixed at 0.4 with G varying. Plots (bl) and 
(b2) have G fixed at 0.01 and e varying. The resulting fidelity after a single successful round, F' , is 
shown in the (al) & (bl) plots. Note that, in the (al) & (bl) plots the LoMM and Zinf protocols 
give identical resulting fidelities, which is shown by a grey curve overlaid with black dots. The 
(a2) & (b2) plots give the overall success probability, P''" = (pfiltcr j2pdistil^ which combines the 
probability of local filtering twice and distilling. 



3.1. From photon loss to generic noise 

Experiments have demonstrated entanglement over macroscopic separations via 
measurements on photon£IEH Building on this success, a number of promising 
schemes for quantum computation propose projecting matter qubits into an en- 
tangled state depending on photon measurement j^HlZI^^I^^P^P^P^P'^P'^PlI However, 
some schemes are not robust against photon loss, which can cause the projected 
state to contain a large proportion of the |1, 1) component. Here we consider mixed 
states damaged by a combination of photon loss noise and depolarizing noise: 

PPL = (1 - e)\^+){^+\ + 1>(1, 1| (24) 

When G = 0, there is only photon loss noise, and the mixed state will locally 
filter to a rank deficient state. When G = 1, all noise is depolarizing and ppl is 
LoMM without any local filtering. For < G < 1 we have the intermediary case 
and the base pair can be locally filtered into a LoMM canonical form. Furthermore, 
the base pair inherently has all local information in the Z-basis, and so Zinf can 
be performed without the filtering step. Distillation of photon loss noise has been 
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relaxation probability: y relaxation probability: y | I I 

Fig. 3. One round performance of distillation of photon loss noise (see Eq. ( |25| ) for various 
protocols. Plots (al) and (a2) have the relaxation rate 7 fixed at 1/2 and vary d. Plots (bl) and 
(b2) have the damping angle d fixed at 7r/8 and 7 varying. The resulting fidelity after a single 
successful round, F' , is shown in the (al) & (bl) plots. The (a2) & (b2) plots give the overall 
success probability, P**" = (pfilt<=r^2pdistil^ which combines the probability of local filtering twice 
and distilling. 



considered previousljEHH and the best-known protocols are instances of our more 
general Zinf protocol. 

In Fig. ([2| we show the performance of four different protocols at distilling 
Ppl- With the noise rate, e, held constant (e = 0.4) and G varying, we see that 
as G — > 1 all protocols converge to the same output fidelity. In this limit, we also 
observe that the overall success probability, P''", converges to the same value for 
all protocols except Zinf, which suffers an inherent 1/2 overhead. However, we are 
really interested in the strong local information regime, G ^ 0, where photon loss 
dominates. In this regime Zinf and LoMM protocols produce exceptionally high 
fidelity states after only a single successful round of distillation, with Horodecki 
distillation and DEJMPS performing significantly worse. Whilst Zinf and LoMM 
give identical outputs on success, LoMM has an overall success probability that 
vanishes with G as Pll^yi ~ (1 - e)^G2 + 0[G% 

Although our protocols achieve higher fidelities in a single round this comes at 
the cost of lower success probability, which can be a significant cost for LoMM. 
In section [4] we consider one approach to applying many rounds of distillation, 
and find that the yield, which depends on both fidelity and success probability, is 
greatest for Zinf in the regime of interest. 
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3.2. Amplitude damping 

When the logical states of a qubit correspond to the excited and ground state of a 
physical system, it can gradually decay down to the ground state. Such a process 
is known as amplitude damping, and despite its prevalence in the laboratory the 
present author is not aware of any studies of distilling entanglement from amplitude 
damped states. Furthermore, since amplitude damping increases local information 
it is a natural choice as a test bed for our proposed protocols. Hence, we consider 
the family of base pairs: 

PAD= (^a<^K^b)\'^+){^+\{K\^K^^)1 (25) 

Here K-^ and Kg are local Kraus operators acting on qubits A and B, such that: 
K - U! (|0)(0U + yr^|l)(lU) U-', (26) 

where 7 is the relaxation probability, and the unitary is: 

[/f = cos(6')l:r + ism{9)Y^. (27) 

When 9^0, the base pair pad describes computational basis amplitude damp- 
ing acting on the l^*^) state. When 6* 7^ the effect of the noise channel can be 
interpreted as either amplitude damping in some non-computational basis, or equiv- 
alently computational basis damping that has acted on some state {U^ ® U^)\'i'~^). 

Figure ([S]) depicts the performance of our protocols and Horodecki distillation 
for a range of parameters, 7 and 9. After a successful first round, both LoMM and 
Zinf always produce states with a greater fidelity than DEJMPS or Horodecki dis- 
tillation. This is most pronounced when amplitude damping approaches the compu- 
tational basis, 9 0, where both our protocols approach an output of unit fidelity. 
In contrast, DEJMPS approaches a separable state in this limii|^as for 9 = and 
7 < 1/2 the initial base pair has a maximally entangled fraction below 1/2. In the 
same limit, we observe the now familiar feature of the LoMM protocol having a 
vanishing success probability, whereas all other protocols maintain a finite success 
probability. Again our protocols offer increased fidelity in exchange for reduced suc- 
cess probability. In section |4] we consider the overall balance of these effects, where 
we will find that Zinf has the greatest yield in the small 9 regime, and all other 
protocols converge to a similar performance as — > 7r/4. 

Distillation of amplitude damping noise also presents some interesting features 
that we did not encounter for photon loss noise. Most notably, the Zinf protocol 

''This is because of our choice of 7 = 0.5. For 7 < 0.5, DEJMPS would produce an entangled state 
for all values of 9. Whereas for 7 > 0.5 there exists a non-zero threshold value of 9, belowr which 
DEJMPS produces a separable state. When the parameters, 7 and 9, are such that DEJMPS 
produces a separable state, the initial base pair has a maximally entangled fraction below 1/2, 
formally Fmax < 1/2. 
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Damping basis angle: 9 Relaxation probability: y 



Fig. 4. The yield (the average number of Bell pairs of fidelity F = 0.99 produced per base pair 
consumed) for various protocols when the base pair suffers (a) photon-loss/depolarizing noise (see 
section |3.1| ), and (b) amplitude damping (see section |3.2| ). Local information is greatest for 
(a) large e and small G, or (b) large 7 and small d. Zinf (dotted line) is the only protocol that 
consistently has the largest yield in these regimes, often gaining over an order of magnitude over 
prior protocols. Note that the ranges of yields covered is not the same across the four figures, and 
in particular (al) and (bl) only go up to a yield of 0.1. 



exhibits a striking discontinuity as 9 varies: it jumps between producing the same 
output as LoMM and producing a less entangled output. This discontinuity occurs 
because we have specified that Zinf uses the local filter that maximizes the quantity 
-^zinf • "^^^ ^^'^ different branches of the discontinuity correspond to local filters from 
different families: (a) the family of local filter that produces pzinf with the variables 
Xi ordered optimally (e.g. Ai largest, A4 smallest); or (b) the family of local filters 
with a suboptimal ordering of variables. However, the reader should note that 
the choice of Zinf local filter is flexible, and for some tasks it may benefit from local 
filters with optimally ordered A^. 



4. Symmetric Scheduling 

So far we have only considered one round of entanglement distillation, which may 
not deliver the required fidehty. Typically we need to repeat many rounds of distil- 
lation. However, the optimal means of repetition depends on a number of factors, 
including: the number of qubits at each location; whether the procedure forms a 
component in a repeater architecture^ and whether memory effects are importanlP. 
The simplest scenario is when there are no size limits on the local Hilbert space, 
memory effects are neglected and we aim to produce Bell pairs between two fixed 
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locations (rather than via several "hops" in a repeater architecture). This simple 
scenario admits a simple, symmetric, approach to repeating distillation. In sym- 
metric scheduling, many copies of the base pairs are simultaneously stored, and 
simultaneously distilled to produce an output denoted p2- Pairs of p2 are then 
distilled into an output pa, and so on. 

For Zin:^ and LoMM, the output of the first round is a Bell diagonal mixed 
state that no longer contains local information. Therefore, our protocols will only 
be used in the first round of symmetric scheduled approaches, with subsequent 
rounds performed by DEJMPS. Nevertheless, Fig. |4] shows a marked difference 
in the yield for all protocols. The yield depends on both the one round fidelity 
and success probability, and so the yield results should be compared alongside 
the one round performance. Specifically, low yields occur whenever the one round 
success probability vanishes or the one round fidelity approaches 0.5 (the value for 
a separable Bell diagonal state). 

We see that for Zinf the reduced success probability can be worth the increased 
fidelity. Indeed, in strong local information regimes Zinf can produce yields 10 times 
greater than prior protocols. Our results show that LoMM can sometimes match 
Zinf 's yield in regimes of modest local information. However, LoMM has a rapidly 
vanishing yield when the local information becomes very significant ( e.g. G ^ or 
6 —*Q). This yield crash of LoMM is due to a rapidly decreasing success probability 
for the local filter into LoMM form. 

In contrast, as 6 7r/4or G — > 1, local information decreases and all pro- 
tocols converge to similar yields. Although yields are similar in this regime, the 
Zinf protocol can perform slightly worse. This is because Zinf postselects on fewer 
measurement outcomes, and so has a factor of 1/2 overhead. 

5. Entanglement Pumping 

Although symmetric scheduling offers a rapid means of producing high fidelity Bell 
pairs, it requires a large Hilbert space to simultaneously store many copies of p. 
This spatial cost can be very large, and in many instances the available technology 
may not be able to provide such a large Hilbert space. In these scenarios, it is ad- 
vantageous to perform distillation between noisy Bell pairs that are not identical. 
Such approaches are asymmetric, as in each distillation attempt one participating 
Bell pair will have a higher fidelity than the other. For example, the most asym- 
metric of approaches is entanglement pumping where one of the two pairs is always 
a base pair. The advantage of entanglement pumping is that we only have to store 
two pairs at any one time. However, the price of conserving Hilbert space is a slower 
convergence to higher fidelities, and a reduced yield. Since entanglement pumping 
is an important application, in subsection |5] we show that Zinf can be modified to 
this task. In our original account of Zinf it was performed between two identical 

^Note that the Zinf filter has been chosen to maximize yield. 
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Bell pairs, and the local information cancelled out. Here we show that the local 
information still cancels out provided that we pump for an even number of rounds. 
By using several levels of entanglement pumping one can achieve higher fidelities. 

In subsection [5] we give a brief overview of the minimal number of pumping levels 
required by prior proposals. We compare this with our pumped-Zinf protocol, which 
proves to be conservative in its spatial requirements. 

5.1. Overview of few-qubit pumping schemes 

Diir and BriegeF* proposed a method of using many pumping levels to enable quan- 
tum computing with only a small amount of spatial resources. In the first pumping 
level, two copies of the base pair p are distilled by a single round of DEJMPS to 
produce p2. Still in the first pumping level, we perform DEJMPS using one copy of 
p and one copy of p2, which upon success produces a mixed state p^. This process 
is continued using one copy of p and one copy of pi , until after n rounds the fidelity 
plateaus for p„. This resource, p„, then plays the role of the base pair in the next 
level of pumping. If we denote p„ — p', then the next level of pumping proceeds 
in the same fashion with p' replacing p. After some m rounds in the 2""^ pumping 
level, we have the state p'„^. This then forms the resource in the 3"^ pumping level, 
such that p" — and so on. 

The Diir and BriegeP proposal uses N + 2 qubits per location to support N 
pumping levels. Of these qubits, A'^ + 1 play an active role in distillation, and 
one qubit stores a logical qubit involved in the computation. The overall strategy 
is to distill a high fidelity Bell pair within the pumping levels and then use it 
to implement a two-qubit unitary between two separated logical qubits. Diir and 
Briegel found that two or three pumping levels, and hence four or five qubits per 
location, attained a high enough fidelity for fault tolerant quantum computing. 

Jiang et al^ proposed a slight variation on pumped-DEJMPS that, assuming 
perfect local operations, enables quantum computing with any entangled LoMM 
base pair using only two pumping levels, and hence four qubits. Base pairs of 
LoMM form with a vanishingly small 4*'^ eigenvalue, such as phase noisy states, 
can be distilled using only a single pumping level. If the available Hilbert space is 
smaller still, then by employing measurement-based quantum computing (MBQC) , 
the same base pairs can enable quantum computation with the use of one less qubilP. 
This one qubit saving is possible because the entanglement resource for MBQC, e.g. 
a cluster state, can be grown by probabilistic parity projectionaSMMZEll 

5.2. Pumped-Zinf 

Whilst the Hilbert space requirements are well documented for LoMM base pairs, 
other mixed states have been neglected. Hence, it is informative to consider the 
Hilbert space requirements of applying LoMM, Zinf or Horodecki distillation to 
non-LoMM mixed states. All these protocols begin by local filtering the base pair. 
However, most physical systems do not admit direct POVM measurements for local 
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Stepl Q P Q Steps 2-5 




(^) Repeat ^ {a)=£ia) 

Steps 2-5 

Fig. 5. A schematic of pumped-Zinf (a similar pumped version of LoMM is also possible). En- 
tanglement pumping uses base pairs, p, to increase entanglement of a state cr; e.g. via noisy parity 
projections, . The state rr may be the output of a previous round of Zinf, e.g. the state p from 
Fig. [l] We always perform pumped-Zinf an even number of times, consuming an even number of 
copies of p. An even number of applications is important because one application of M will pass 
local information from p onto cr, but with two applications of Ai this local information cancels 
out. Since pumped-Zinf creates an imperfect parity projection, it is possible to entangle qubits A 
and B of cr without destroying prior entanglement held with third parties. We illustrate that cr 
may hold third party entanglement by lines emerging from cr, and this has a natural interpretation 
in terms of graph state growth. Unlike in Fig. [T] we do not begin by preparing a locally filtered 
Pzinfi and hence we eliminate the need for an additional ancillary qubit that would be required 
to perform the POVM measurement. 



filtering, and so they have to be implemented indirectly. Fortunately, it is well known 
that a single ancillary qubit can be used to implement any two-outcome POVM 
measurement. Hence, one may expect that distillation of general mixed states may 
require at least one more qubit per location than required for distilling LoMM base 
pairs. 

However, this extra qubit is not essential. Instead of using an ancilla to attempt 
the local filter, the same resulting quantum operation can be achieved by a modified 
protocol. Our modified protocol uses only local unitaries and computational basis 
measurements and so requires no greater technological sophistication than prior 
protocols. For completeness, we will also modify our protocol to distill a high fidelity 
parity projection, so it can be used to either produce Bell pairs or grow entangled 
resources for MBQC. We describe an entanglement pumping version of Zinf where 
each round involves a copy of p, and some other state a. As before p can be related 
to a state, pzinf, by local filters gA (8> gs- Later we will use that these local filters 
can be expressed as 

gA^ao\+){Ao\A + ai\-){Ai\A, (28) 

93 = Po\+){Bo\B + fh\-){Bi\B, 

where \Ai) are normalized, but not necessarily orthogonal, and similarly for \Bi). 
The variables and Pi are complex numbers satisfying ja^l, |/3i| < 1. Our pumped- 
Zinf protocol requires no additional ancilla qubits to perform local filtering, and 
is capable of distilling probabilistic parity projections for graph state growth (see 
Fig. [5] for outline). We use Al and Bl to denote qubits holding one copy of p, 
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and A2 and B2 as the qubits in the state a that we wish to parity project. Our 
pumped-Zinf protocol is then: 

(1) Prepare qubits Al and Bl in the base pair state p; 

(2) Perform local unitaries Ua®Ub \ 

(3) Measure qubits A\ and B\ in the computational basis, post-selecting on the 
outcome |1)ai|1>bi; 

(4) Perform local unitaries Va®Vb \ 

(5) Measure qubits A\ and _B1 in the computational basis, post-selecting on the 
outcome 

(6) Repeat for subsequent rounds an even number of times; 

The protocol is quite complex and involves a number of new unitaries: 

Ua = (|l)(Ao| + |0)(lo|)^^ ® |0)(0U2 + + |0)(li|)^^ ® |1)(1U2, (29) 

Ub = (|l)(So| + |0)(So|)^^ ® |0)(0|b2 + + |0)(Si|)^^ ® |1)(1|b2, 

Va = W{ao)Ai ® |1)(0U2 + W{ai)Ai ® |0)(lU2, 
Vb = W{P^)bi ® |1>(0|b2 + W{Pi)bi ® |0)(1|b2, 

where W{c) can be any one qubit unitary satisfying (l|VF(c)|l) = c. The barred kets, 
e.g. I^o), are orthogonal to their unbarred counterpart so as to ensure unitarity. The 
role of these unitaries in pumped-Zinf can appear opaque at first, and so we pause 
to relate them to something more familiar. Consider local operations at location A 
for steps 2-5 and we find: 

{1\aiVa\1){1\aiUa = aoMoUi ® |1)(0U2 + «i(AiUi ® |0)(1U2, (30) 
= XA2{'^\CX^lgAi, 

and so these local operations emulate the behavior of performing the desired local 
filter, followed by a control-not and post-selected measurement. Hence, this involved 
protocol manages to reproduce the effect of local filtering and distillation without 
the use of an additional ancillary qubit! Actually, a local filter is performed, via 
the POVM measurement in steps 4 — 5, which reuses qubits A\ and Bl after they 
where measured in step 3. Our protocol can be thought of as a temporal reordering 
where local filtering is performed after distillation, rather than before, hence saving 
valuable spatial resources. 

A single successful round of the protocol generates a quantum operation, which 
we denote M{a). In each application of A^(ct) the presence of local information 
has the effect of partially projecting a into some separable state. However, the 
preferred subspace changes for alternate rounds, cancelling itself out after an even 
number of applications of M.. This cancelling of local information emerges from 



continual bit-flipping, which is explicitly presented on the R.H.S of Eq. (30). If 
these bit-flips where not present then local information would accumulate. Hence, 
we always apply our protocol for an even number of rounds, and actually work with 
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£{<t) = M{M{(t)), where: 

£{a) « fi+P+ ((1 - q+)a + q+ZA2^ZA2) P+ (31) 
+^1^P- ((1 - q-)a + q^ZA2<yZA2) P", 

where P+ {P^) are even (odd) parity projectors, P"^ = (1 ± Z^2^B2)/2. As with 
the original Zinf protocol, this is related to the concurrence eigenvalues or p, such 
that: 

/i+ = (Ai + A4)' , q+ = 2AiA4/Ai+, (32) 
Ai- = (A2 + X^f , q- = 2X2X3/ fi-. 

Expressing the quantum operation in terms of these variables proves useful when 
calculating the effect of repeated pumping. It is helpful to compare these results 
with the analogous Eq. ( 23 1 , where the terms match ujF] We can also calculate the 



Jamiolkowski fidelity^Sl of £ compared to a perfect parity projection P+ , and one 
finds that one application of £ has the same fidelity as one round of the unmodified 
Zinf protocol (see Eq. (14)). 

Now consider 2n successful rounds of pumped-Zinf, which is equivalent to n 
applications of the quantum operation £: 

£"{a) cx ^ilP+Vl{a)P+ + ^ilP-V'l{a)p- , (33) 

where I'5((t) represents n applications of a dephasing channel: 

T^±{'^) ^ - q±)<^ + q±ZA2crZA2- (34) 

Examining the quantum channel f "(cr) we can see two separate effects at work: (a) 
as n increases the odd parity projection component will vanish (assuming > 
/i_ ) ; (b) the quantum operation becomes more dephasing. In fact the Jamiolkowski 
fidelity of factorizes with respect to these separate effects: 

-^(t(:)<-..r-'ri). (-) 

Such behavior is not unique to our protocol but occurs in other pumped protocol^3H 
Earlier we commented that some protocols only require a single pumping level to 
distill certain Bell diagonal states to unit fidelity. These special, easily distillable, 
states need to have only one vanishingly small eigenvalue. Here we have a much 
stronger result: if a base pair p has a vanishingly small concurrence eigenvalue 
A4 = 0, then it can be distilled to unit fidelity with only a single pumping level. 
We can distill towards unit fidelity because there is no dephasing in the even parity 
subspace, since A4 = entails q^ = 0, which gives the simple fidelity expression: 

^(^4-0) = -^. (36) 



^To see this one has to multiply out the new variables, such as ^t+(l — g+) = + A|, to see they 
equate to the earlier expressions. 
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This result substantially broadens the class of mixed states that can enable quantum 
computing with only two qubits per location. 

What if a base pair does not have a small enough concurrence eigenvalue? In 
this instance it seems unavoidable that an extra pumping level is required to enable 
high fidelity distillation. When pumping Zinf with two levels, we prepare a in the 
state I +,+)(+,+ 1 so that ^"((t) is a LoMM mixed state. This LoMM mixed state 
is then used in the next level of pumping with no local filtering required. 

6. Conclusions 

We have proposed two protocols, LoMM and Zinf, for entanglement distillation. 
Our protocols can achieve high fidelities in fewer rounds of distillation than previ- 
ous protocols. The mixed states we produce have fidelities that are a function of the 
concurrence eigenvalues encountered when calculating the concurrence of the base 
pair. Our protocols achieve this by exploitation of local information, which is dis- 
carded in DEJMPS. Although LoMM can attain high fidelities in a single round, in 
the strong local information regime it suffers from a vanishing success probability, 
and hence yield. Whereas Zinf can achieve the same fidelity for only a small over- 
head in success probability. Consequentially, symmetrically scheduled Zinf achieves 
much higher yields for base pairs with large amounts of local information. 

We have also shown that Zinf can be modified into a compact entanglement 
pumping protocol. Furthermore, pumped- Zinf enables measurement-based quantum 
computing in distributed architectures with a very small number of qubits per 
local site, for any base pair. In general pumped-Zinf requires 3 qubits per local 
site. However, only two qubits per location are needed when the base pair has a 
vanishing concurrence eigenvalue. 

Zinf offers increased yields when spatial resources are abundant and can function 
when spatial resources are very constrained. Since Zinf has much to offer on both 
extremes of the spatial scale, we also expect that it will have applications in the 
intermediate regimes. 
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sidered. This work has been funded by the Royal Commission for the Exhibition 
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